Physical Review Research

The relationship between genotypic and phenotypic variation is determined by the complex interaction of genetic and environmental factors. While statistical methods capable of detecting such interactions exist, an axiomatic mathematical framework that seamlessly describes the combined effects of genetic modifications and environmental exposures on a common scale is lacking. In this report, buffering concepts are used to construct a measurement system that enables the geometric representation of both gene-by-gene and gene-by-environment interactions on the extended complex plane (i.e., as projections on the Riemann sphere). In this manner, any such interaction, or combination thereof, can be precisely defined and quantified as the deviation from the neutral value calculated through the applicable complex transformation. When thus conceptualized, the framework's parameterization defines the "state space" of a given measurable phenotype along both the real and imaginary dimensions, thus establishing an unambiguous and broadly applicable method for determining the phenotypic value expected upon combinatorial changes in genetic and/or environmental variables. Remarkably, by applying these methods, it is possible to quantify the effects of any gene-by-environment interaction using the equation, AGxE=Im([z]obs*zexp)/2, where zobs and zexp are complex numbers representing the observed and expected phenotypes of a given genotype expressed in terms of the buffering parameters, and b.

Diffuse gliomas are deadly because the individual tumor cells invade - they travel far from the imageable mass, so it is impossible to remove the tumor completely. On the cellular level, glioma cells seem to be in either a "go" state (in which they do not divide) or a "grow" state (in which they do not migrate). We investigate what this tiny choice has to say about the large-scale speed of the invasion front and whether the implication is sufficiently strong to rule out the classical description of the Fisher-Kolmogorov-Petrovsky-Piskunov (Fisher-KPP) type, in which a single phenotype migrates and proliferates. We derive a two-phenotype reaction-diffusion model with density-dependent switching, and we prove the cooperative (quasi-monotone) structure and the associated comparison principle and study travelling-wave solutions of the model. A leading-edge linearization gives minimal front speed as minimizer of an explicit dispersion relation, and direct simulation verifies the predicted speed. In the experimentally relevant fast switching limit, we find a closed-form expression for the speed, that is, we obtain an effective Fisher-KPP equation with rescaled diffusivity and growth rate, with the fractions of the phenotypes. The "go-or-grow" (GoG) front can move at a maximum speed of half the Fisher speed for the same single-cell motility $D$ and proliferation rate $r$, which occurs only when the cells divide their time equally between the two phenotypes. This bound is directly testable: measurement of the front speed, plus independent determination of $D$ and $r$, discriminates the two hypotheses, and in the GoG case, yields recovery of the phenotype balance. We then extend the result to anisotropic (DTI-informed) invasion along white-matter tracts and discuss implications for understanding clinical measurements of growth rate.

Biomolecular condensates are widespread cellular self-assembled structures with essential functions. There are suggestions of condensates formed by different proteins being near criticality. However, systematic investigation of the criticality of condensates is absent, and critical exponents defining their universality class have not been found. Here, using long-time simulations, we show that condensates exhibit typical critical phenomena, including scale-free spatiotemporal correlations, critical slowing down, divergence of correlation length and dynamic scaling. From these scaling behaviors, a set of critical exponents is determined. Based on dynamic critical exponent, diverse condensates can be divided into two distinct universality classes, arising from differences in their molecular components and interaction types.

Fixational eye movements (FEMs) comprise the involuntary small scale eye motion conducted during fixation on a stationary stimulus. As a consequence, the visual information can be spread across multiple photoreceptors reducing the local signal-to-noise ratio. Yet, the signals transmitted by individual photoreceptors adapt to constant stimulation so that an entirely still scene would eventually fade from view. Because FEMs convert a stationary stimulus in the world to a temporally varying one on the retina, they can act to prevent this stimulus fading. Thus, FEMs can be understood as a sampling protocol than needs to be adjusted to the underlying processing circuitry. We analyse the impact of FEMs on the rate of information acquisition at the level of the retina for two common tasks of the human eye that typically go hand in hand: detection and localization. Here, we build a simple analytical model of visual perception, i.e. we subject a continuous receptor array to a stimulus moving across the retina as a consequence of FEMs with receptor excitations depending on past stimulation through a linear response function. Using Bayesian inference we quantify both the probability of detection and the accuracy of localization as a function of parameters controlling eye movements and stimulus. We find that localization of a stimulus is equivalent to the detection of the stimulus gradient. This allows us to discern optimal properties of eye movements for the respective tasks and provides a link between two typical psychophysical observables: detection thresholds and Vernier acuity. Our analysis suggests that typical human FEMs tend to facilitate localization at the expense of detection. Simply put, if you can see a stimulus you also know where it is. Finally, we propose a variety of experimental protocols to investigate the interplay between FEMs, detection, and localization with the potential of inferring intrinsic properties of an individuals visual system.

Mechanical properties of epithelial tissues play essential roles in morphogenesis and physiological function. In this study, we analytically derived the in-plane bulk modulus, shear modulus, and Poissons ratio of a three-dimensional cell vertex model of epithelial monolayers. We showed that the model can robustly reproduce a near-zero in-plane Poissons ratio, a mechanical feature reported in cultured epithelial tissues. Numerical simulations further confirmed that the theoretically predicted Poissons ratio accurately describes the response of the model under finite, biologically relevant strains. In addition, the model exhibits not only morphological bistability between squamous-like and columnar-like states, but also mechanical bistability characterized by distinct elastic responses. Together, these results provide a minimal three-dimensional framework that links cell-scale mechanical interactions and epithelial morphology to tissue-scale elastic properties.

Since the earliest microscopic observations, the geometric organisation of cells has captured biologists interest. Recent work by Gorgi et al. showed that bacterial colony organisation, including biofilms, can be explained across diverse species by radial expansion from fixed initial seeding sites and contact-inhibited growth, with little need for species-specific mechanisms. Here, we extend this geometric framework by incorporating seeding time as an additional driver of colony organisation. Using simulations and analytical models for expected colony size, we show that staggered seeding yields order of magnitude increases in the expected size of early seeded founder colonies. At realistic biofilm growth rates, a 2-day lag between founder and subsequent colony seeding produces an approximately 10-fold increase in expected founder size, while a 1-week lag produces a 25-fold increase. These findings provide a simple geometric basis for biological priority effects, illustrating temporal advantage alone can generate substantial spatial dominance, with implications for cardiovascular devices where host and bacterial cells compete in a race for the surface.

We develop a theoretical proposal linking vacuum stability and brain dynamics through superconductivity-inspired coherence, symmetry reduction, and the thermodynamic stabilization of low-entropy regimes. We take an unbroken SU(3) structure as a candidate stable residue of the low-temperature vacuum. At the neural level, we formulate a coarse-grained analog in which a two-fluid model with dissipative and coherence-supporting components describes brain dynamics. Specifically, the coherence-supporting component is proposed as a possible basis for the efficient binding and integration required to sustain a stable, unified conscious state. The proposal offers a common geometric language for relating physics and neuroscience with falsifiable signatures in coherence and state-dependent transitions. The main technical contribution is a computational algebraic model of conscious-state dynamics, where neural data are mapped to reconstructed state trajectories. Effective generators are inferred from those trajectories, and the two-fluid split is tested as a Cartan-root decomposition of su(3), with a rank-two commuting sector for coherence-preserving balance and six root directions for state transitions. This structure can be tested on neural data and contrasted with alternative dynamical models.

Touchscreens have become the dominant interface in consumer electronics, yet interactions with them remain primarily visual. Incorporating haptic feedback that simulates touch sensations could make these interactions more natural and intuitive. Electrostatic actuation, which modulates friction by attracting the finger toward a capacitive surface using an alternating voltage, offers a promising approach. The resulting increase in friction is often attributed to the rise in real contact area; however, direct experimental evidence linking voltage input parameters to real contact area and contact forces remains limited. Here, we use frustrated total internal reflection to directly image the real contact area while simultaneously measuring contact forces during controlled finger sliding under electrostatic actuation. We systematically vary voltage amplitude (75-150 V) and excitation frequency (30-230 Hz) and quantify the changes in contact area and forces as functions of these parameters. Our results reveal a quadratic dependence of real contact area, electro-static attraction, and tangential force on voltage amplitude, with comparatively small effects of excitation frequency. These findings clarify the respective roles of voltage amplitude and frequency in the electrostatic modulation of finger contact mechanics, providing design guidelines for haptic display design.

Reentrant phase transitions, in which multicomponent systems phase separate at intermediate concentrations but dissolve at experimentally accessible higher concentrations, are ubiquitous in mixtures such as biomolecular condensates. We show that introducing reversible dimerization into a multicomponent Flory-Huggins model and integrating out the dimer state generate an effective three-body repulsion that reshapes phase diagram geometry. This emergent higher-order interaction arises naturally from an interaction competition and mass-action equilibrium, providing a microscopic explanation to the previously phenomenological three-body interaction. We derive a closed-form phase boundary equation capturing phase separation and reentrant dissolution in this minimal model, and predict explicit interdependence between competition strength, emergent many-body interactions, and dissociation constants. We recover and extend the reentrant phase boundary scaling relations through interaction renormalization, with regime of validity. We apply our model to G3BP1-RNA-suramin and explain the underlying mechanisms from the physical parameters inferred from reentrant phase boundaries.

Transcription factors organize into liquid-like condensates to facilitate gene expression, yet the physical mechanisms governing their formation and properties remain poorly understood. We study the size statistics of transcriptional condensates in human HAP1 cells using widefield and super-resolution microscopy tagging the epigenetic reader BRD4. We find that hubs that appear monolithic in widefield resolve into clusters of smaller droplets that resist coarsening. We link this size control to Active Model B+, a non-equilibrium field theory that captures a regime of reverse Ostwald ripening out of thermal equilibrium. In this regime, chemically driven currents cause larger droplets to transfer mass back to smaller ones, stabilizing a state of microphase segregation. The observed exponential size distribution of BRD4 foci quantitatively matches our numerical simulations, suggesting a universal physical picture for the non-equilibrium self-limitation of cellular condensates.

Brain dynamics can be described in three different convenient mathematical languages, namely connectome harmonics, turbulence and complex harmonics (CHARM). Here we demonstrate that these theoretical frameworks can be rigorously unified, under the functional calculus, as one self-adjoint operator and its single spectral measure. The connectome Laplacian carries that measure; the harmonics are its spectral projections, the turbulence smoothing kernel is its resolvent, and the CHARM form is its unitary propagator. The bridge that makes this exact is a textbook fact: The exponential distance rule, which is the empirical kernel of the turbulence model, is the Greens function of a screened Laplacian, so the local order parameter is the phase field passed through the resolvent of the same operator whose eigenfunctions are the harmonics. A single shared control parameter, the spectral gap, simultaneously yields the cortical hierarchy, the turbulent information cascade and the structured interference the CHARM form measures. This unification makes a strong predictive claim. If the harmonic projections, the turbulence resolvent and the CHARM propagator really are three functions of one operator, then any structural perturbation that re-tunes the operator must move all three signatures in unison and must do so with a single coupling. We test this prediction with a pharmacological perturbation by lysergic acid diethylamide (LSD), which is known to change the emotional state, by empirically perturbing the operator with a 5-HT2A receptor density map and asking whether one scalar coupling can simultaneously predict the multi-scale turbulence shift observed, through the resolvent, and the macroscale harmonic energy redistribution, through first-order Rayleigh-Schrodinger perturbation theory. We found that the two independent functional domains respond in unison to one structural perturbation of one operator. The identity is exact as operator calculus and its purchase on the brain depends on a single load-bearing seam, the degree heterogeneity of the connectome, which we make explicit. We propose that this single-operator structure is the necessary mathematical scaffolding of our Entangled Loop theory.

Cell migration is fundamental to various biological processes, including morphogenesis, wound healing, and cancer metastasis. Durotaxis--directed migration of cells in response to spatial variations in stiffness--has been extensively studied using engineered substrates with prescribed stiffness. However, recent work has increasingly shifted toward understanding cell migration in fibrous matrices that can be actively remodeled by the actomyosin contractility, as commonly observed in tumor and epithelial cells. Despite these advances, a theoretical framework explaining how cells structurally remodel their surrounding matrix to promote their own durotaxis, and which cellular forces govern this behavior, remains elusive. To address this gap, we developed a biomechanical model in which polarized cells contract and migrate over a fibrous matrix. Using this model, we first confirmed that cells on an externally strained matrix preferentially migrate along the direction of applied strain. Then, we investigated how cells autonomously remodel the matrix to create stiffness patterns favorable for durotaxis. In the presence of intercellular adhesion, cells acted collectively to stiffen the matrix, after which a small subset of cells escaped the main population and migrated outward. This behavior is reminiscent of intravasation during cancer metastasis, where cohesive cell clusters generate local matrix remodeling that facilitates the departure of more motile subpopulations. These results illustrate how matrix stiffening driven by cell cohesion and contractility regulates durotactic behavior and provide mechanistic insight into collective invasion processes relevant to cancer metastasis.

Cetaceans reduce the net cost of sustained swimming through intermittent locomotion, alternating active fluking with unpowered gliding. The energy balance of this strategy is central to understanding survival rates, population sustainability, and the effects of anthropogenic and environmental pressures. While active-phase energetics have been characterized extensively, the glide phase remains largely unexplored. Here we derive the optimal glide duration (Topt) and the maximum glide duration beyond which energy savings vanish (Tzero) for three odontocetes spanning a 20-fold range in body mass, using high-fidelity CAD models and wall-modeled large eddy simulations. We show analytically that speed retention at Topt and mass-specific peak energy savings are both fully determined by the active-to-passive drag ratio, propulsive efficiency, and swimming speed, independently of body morphometry and drag coefficient, and are therefore invariant across species at any given speed. These passive-phase optima extend the known size-independent active-phase invariants to the glide phase, towards a scale-independent energetic framework for burst-and-glide locomotion in small cetaceans.

Understanding how complex systems self-organize, exhibit emergent properties beyond their constituent elements remains a challenge across physics, biology, and cognitive science. In resource-constrained neuronal systems, existing theoretical approaches, including gauge theoretic formulations, statistical physics-inspired methods, dynamical population models, and variational principles such as the Free Energy Principle, address important aspects of this problem but do not fully specify the physical conditions and thermodynamic costs under which self-organizing behavior occurs. Here, we introduce Dynamic Resource Theory (DRT) as a general physical framework for describing self-organization under constrained resource availability. DRT formalizes complexity as a physical property of self-organizing systems arising from coupled mechanisms of resource allocation and dynamic reallocation of internal resources. This framework provides a thermodynamic and variational account of how stability is preserved while adaptive reconfiguration remains possible, consistent with stationary action and thermodynamic constraints. DRT is formulated within a gauge theoretic setting and directly incorporates the energetic costs associated with maintaining structure and enabling system-level reconfiguration. Within DRT, baseline resource allocation preserves system stability, while internal and external demands perturb the system, driving self-organization through dynamic resource reallocation across a coupled free energy landscape without assuming subsystem separability. We then develop Neural Resource Theory (NRT) and Cognitive Resource Theory (CRT) as principled specializations of DRT, illustrating how this structure is instantiated in resource constrained neuronal and cognitive systems. We conclude by discussing the broader implications of DRT for understanding how complexity, emergence, and adaptive capacity arise over time through thermodynamically permissible reallocation processes across scales.

The hierarchical organization of multiphase biomolecular condensates into core-shell architectures is a fundamental problem in soft matter and biophysics. While classical explanations rely on hierarchies of interfacial tension ({gamma}) between coexisting liquids, the ultralow tensions of condensates (0.1-1 {micro}N/m) render such hierarchies potentially fragile. We introduce a robust assembly principle based on Polymer-Assisted Condensation (PAC), in which a single polymer species dictates the entire structure. The polymer nucleates a dense core by recruiting a condensation-incompetent protein (P1). A second incompetent protein (P2), which is repelled or otherwise thermodynamically disfavored from entering the polymer-rich core, is nonetheless recruited to the interface by weak attraction to P1, forming a stable shell. This effective repulsion-driven layering operates across a wide parameter space without requiring{gamma} asymmetries and yields a robust structure that is impervious to concentration fluctuations and environmental perturbations. Phase-field modeling and molecular simulations establish this mechanism and capture key features of nucleolar organization. Our work reveals a general physical pathway for encoding spatial order in soft, multicomponent fluids.

The emergence of amino acids (AAs) and nucleobases (NBs) across meteorites, interstellar ices, and laboratory shock experiments presents a paradox: why do these specific molecular motifs--a minuscule subset of organic chemistrys combinatorial space--appear repeatedly across diverse environments, in the absence of biological selection? We identify a physical mechanism, prebiotic selection, which biases driven chemical systems toward configurations with high stationary probability p*(x) under sustained entropy flux. The bias is quantified by an information quasi-potential {Phi}I (x) = - ln p*(x), entering the overdamped Langevin dynamics O_FD O_INLINEFIG[Formula 1]C_INLINEFIGM_FD(1)C_FD where {Sigma} is the local entropy production rate (Schnakenberg 1976). {Phi}I is defined self-consistently via the full non-equilibrium stationary density, avoiding the circularity of identifying it with a scalar potential. Two central theorems underlie the framework. Theorem 1 establishes that {nabla}{Sigma} and {nabla}{Phi}I are generically linearly independent off equilibrium, so the dynamics is genuinely two-field. Theorem 2 (structural constraints on single-field gradient dynamics) shows that single-field models on compact manifolds (i) produce yield curves that are at most unimodal under linear driving, and (ii) combine disjoint perturbations additively, giving superlinearity factor S = 1 + O(||{delta} V ||2). The observed superlinear synergy of Ferris et al. (1996) lies far outside this perturbative bound and therefore requires the two-field structure of EOM-IFF; the non-monotonic peak of Blank et al. (2001) is consistent with two-field dynamics and also with single-field dynamics in the unimodal-with-peak case of Theorem 2 part 1, so it does not by itself discriminate. From these results, we: (i) define a formal substrate-minimal criterion for prebiotic selection; (ii) show consistency with the non-monotonic shock-synthesis yield of Blank et al. (2001) (R2 = 0.885, peak at P* = 28.4 {+/-} 1.4 GPa); (iii) show consistency with the superlinear clay-catalysed RNA polymerisation of Ferris et al. (1996) (synergy factor S {approx} 5.75, robust under {+/-}1-nucleotide measurement uncertainty); and (iv) state two further falsifiable predictions awaiting dedicated experimental tests. Every lemma and theorem is accompanied by explicit assumptions, regime of validity, and regime of failure; the frameworks scope is what it claims, not more. Prebiotic selection is identified as a physical process distinct from and prior to biological selection, offering a unified account of chemical convergence in carbon-nitrogen chemistry under sustained entropy flux.

Interactions between distinct populations of excitatory (E) and inhibitory (I) neurons can produce complex dynamical landscapes, featuring multistability, oscillations, and paradoxical perturbation responses. By employing an elementary model, the threshold-linear network (TLN), we indicate mathematical conditions for each dynamical regime across fundamental microcircuit architectures, thereby mapping previously unrelated systems neuroscience hypotheses to a common reference space and obtaining novel insights on inputs and connectivity. Namely, we compare balancing strategies in inhibition-stabilized E-I networks, we interpret experiments on gamma oscillations in a canonical neocortical E-I-I circuit, and we discuss bistability in hippocampal E-I-I networks. Then, we show that connectivity determines three fundamentally different kinds of interactions between assemblies in E-E-I circuits. Moreover in, E-E-I-I circuits we find that balanced clustering hinders lateral inhibition, while opponent clustering can produce different bistable configurations, even between completely unstructured assemblies. We conclude that TLNs allow to grasp deep and universal aspects of microcircuit dynamics.

Understanding how network connectivity shapes neural representations is central to systems neuroscience. While dimensionality reduction methods uncover low-dimensional manifold structure in population recordings, a rigorous framework connecting manifold geometry to network mechanisms and information encoding remains lacking. We develop a differential geometric approach for analyzing neural manifolds in rate-based recurrent networks receiving tuned feedforward inputs. We derive expressions for the pullback metric of neural manifolds, showing how input tuning curves, feedforward and recurrent synaptic connectivity shape manifold geometry. Critically, we establish that the Fisher information matrix at steady states also has the structure of a pullback metric, directly linking intrinsic manifold geometry to stimulus discriminability and information encoding. For noise with slow temporal correlations propagated through the network, we show that recurrent effects on information geometry cancel: Fisher information depends only on the feedforward connectivity. Thus, feedforward connectivity critically determines representational geometry. As an example, we demonstrate that the representation of space by a module of hexagonal grid cells is approximately isometric for random distribution of grid phases. Moreover, a linear feedforward transformation can map spatially random input tuning curves into a population of hexagonal grid cells, forming a toroidal manifold. Thus, feedforward connectivity alone can generate structured spatial representations without requiring carefully tuned recurrent connectivity or continuous attractor dynamics. Recurrent connectivity, however, is shown to improve stimulus encoding under fast noise, thereby implementing a selective noise reduction.

Two-dimensional accounts of consciousness that distinguish global integration from functional diversity are empirically supported [1,2] but lack a formal phase-structure: they do not specify the nature of the transition between the two regimes, the order parameter that governs it, or the quantitative predictions that follow. We provide this structure. We propose that the dimensionality of the neural correlation structure, operationalised as the Participation Ratio of the covariance eigenspectrum, constitutes a second, independent order parameter D that governs a phase transition distinct from global integration {Phi}. Formalised through a Landau-Ginzburg free energy functional F[{Phi}, D], this transition defines a Redundant Integrated State {Delta} (high {Phi}, low D) in which globally integrated mental function is present but phenomenal experience is absent, a thermodynamic phase, not a point on a continuum. The framework generates three falsifiable predictions absent from prior work: (i) a power-law scaling D* [~] |{Phi} - {Phi}_c|^{nu} with measurable critical exponent{nu} ; (ii) a diverging susceptibility {chi}_D = {partial}D/{partial}{Phi} at the consciousness threshold, quantifiable from perturbational EEG; (iii) an explicit dissociation between MCS and VS patients in the ({Phi}, D) space, with MCS predicted to occupy state {Delta}. These predictions are directly testable with existing methodology and are not generated by any current theory of consciousness.

Experimental studies show that tumor cells adopt migratory or proliferative phenotypes depending on the local extracellular matrix (ECM). In this work, we propose a minimal go-or-grow invasion model, comprising two specialist cell phenotypes: proliferating and migratory, with phenotypic switching and cell migration depending on local ECM density. Numerical simulations of this model reveal that the spatial arrangement of proliferative and migratory cells depends on the choice of phenotypic switching function. We then ask whether this specialist cell-phenotype model can be reduced to a generalist cell-phenotype model. We derive a relationship between the reduced model and go-or-grow model in the fast phenotypic switching regime. We observe that the reduced model captures the dynamics of the original model, for a range of realistic phenotypic switching functions. We analytically derive the minimum traveling wave speed of the reduced model in a homogeneous ECM bed. Moreover, using linear stability analysis on the go-or-grow model, we recover the same wave speed expression. In addition, we numerically explore how the key parameters influence the traveling wave speed profile. Our analysis indicated the counter-intuitive result that the wave speed is independent of the matrix degradation rate, and our simulations show that, at most, the speed is weakly dependent on this parameter.